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Abstract 

The vortex-vortex and vortex-antivortex correlation functions are determined 
for the two-dimensional 0(2) model undergoing phase ordering. We find rea- 
sonably good agreement with simulation results for the vortex-vortex correla- 
tion function where there is a short-scaled distance depletion zone due to the 
repulsion of like-signed vortices. The vortex-antivortex correlation function 
agrees well with simulation results for intermediate and long-scaled distances. 
At short-scaled distances the simulations show a depletion zone not seen in 
the theory. 

PACS numbers: 05.70.Ln, 64.60. Cn, 64.60.My, 64.75.+g 
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I. INTRODUCTION 



The late-stage ordering of systems in the process of breaking a continuous symmetry is 
dominated by the dynamics of topological defects. In the case of the n-vector model with 
the number of components n of the order parameter equal to the spatial dimensionality d 
one has point defects which are vortices for n = 2 and monopoles for n = 3 [|T|. We focus 
mainly on the case n = d = 2 here. Because of the conservation of topological charge, 
the ordering in these systems is through the charge conserving process of vortex-antivortex 
annihilation. The statistical description of this annihilation process will be a central focus of 
this paper. We present here the first calculation for the separate vortex-vortex and vortex- 
antivortex correlation functions. Knowledge of these functions is an important ingredient 
in understanding the vortex annihilation process. Since these functions contain detailed 
information about the vortex correlations, comparison of our results with simulation and 
experiment provides the most stringent test of the theory to date. In the original work p[ in 
this area only the signed defect correlation function was computed. The unsigned quantity 
is technically much more difficult to evaluate as we discuss in this paper. The comparison 
of the result for the signed quantity with simulations was confused by the use of a theory 
for describing order- parameter fluctuations which led to unphysical singularities at short 
distances in this quantity. The source of this singularity was recently uncovered by us |^ 
and the theory reorganized so as to give physical results for the signed defect correlation 
function. This development motivated us to make a renewed effort to evaluate the unsigned 
defect correlation function, with the results presented here. 

We show in Figure |l| the vortex- vortex correlation function for the two-dimensional 0(2) 
model for simulations and the theory developed here. Beautiful experiments on two- 
dimensionally aligned nematic liquid crystals, which mimic the 0(2) model, have also been 
able to measure the vortex correlations. Their results essentially track the simulation results, 
but are more noisy, so we do not show them here. There is reasonable agreement between 
the theory and the simulation results. Both show a depletion zone at short-scaled distances 
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for like-signed defects. This is expected on physical grounds since like-signed defects repel 
one another. It is non-trivial that we see evidence for a repulsive interaction. The vortex 
interactions emerge naturally in the theory from the equation of motion for the order- 
parameter. In Figure ^ we show the vortex-antivortex correlation function for theory and 
simulation |^. There is good agreement at large and intermediate scaled distances. There 
is a clear discrepancy between theory and simulation results at short-scaled distances. The 
theory shows a monotonic behavior as the separation distance goes to zero. The simulation, 
however, shows a maximum at short separation distances and then falls rapidly to zero. The 
depletion zone seen in the simulation data for the vortex-antivortex correlation function is 
harder to understand physically since the pair is attractive and headed toward annihilation. 
Some possible explanations will be presented in Section VIII. While the theory satisfies the 
sum rule implied by topological charge conservation, it does not appear that this general 
constraint is satisfied by the simulations. 



II. MODEL 

We consider the n- vector model with 0{n) symmetry, which describes the dynamics of a 
non-conserved, n-component order-parameter field ipi'^) = ("^1(1)5 ' ' ' y4'n{^))- Here we use 
the shorthand notation 1 = (ri,ti). The order-parameter evolves via the time-dependent 
Ginzburg-Landau equation 

f ^VV-^ (2.1) 

dt dxfj 

which can be derived from a free-energy containing a square-gradient term and a potential 
term, V[iIj]- We assume that the quench is to zero temperature where the usual noise 
term on the right-hand side of ( |2.1| ) is zero |0. The potential V[iIj] is chosen to have 0{n) 
symmetry with a degenerate set of equilibrium minima at ip = = ipo- Since only these 
properties of Vlip] will be important in what follows we need not be more specific in our 
choice for V[iIj]. It is believed that our final results are independent of the exact nature of 
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the initial state, provided it is a disordered state. Indeed, it is well-established that for 
late times t following a quench from the disordered to the ordered phase the dynamics obey 
scaling and the system can be described in terms of a single growing length L{t), which is 
characteristic of the spacing between defects. In this scaling regime the order-parameter 
correlation function 

C{12)^{^{l)■^f{2)) (2.2) 

has an equal-time scaling form 

C(r,t)=^2^(x), (2.3) 

where is a universal function oi x = r / L{t) (r = |r| = |r2 — ri|), depending only on n and 
d. It is also well-established that, in the scaling regime. Lit) ~ t't'. For the non-conserved 
models considered here the exponent (j) = 1/2 0. 

III. GENERAL FEATURES OF POINT DEFECTS 

Let us briefly review the theoretical picture associated with topological point defects. The 
results of this section are quite general and hold independent of the particular dynamics, 
such as ( ^?T|) , or the approximation scheme, such as the gaussian approximation, used to 
model the system. The vortex charge density can be written in the form 

p(l)=^g, 5(ri-x,(ti)) (3.1) 

a 

where ^aiti) is the position at time ti of the a*'* point defect, which has a topological 
charge q^. We will restrict the analysis here to the case of charge ±1 vortices where = 1. 
This case dominates the late-stage ordering since higher-charged defects are energetically 
unfavourable and dissociate into charge ±1 defects early on. 

The next step is to note, as pointed out by Halperin [^, that the positions of defects 
are located by the zeros of the order-parameter field ip. Therefore the charged or signed 
density for point defects is given by 
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p(l) = 5m)]V{l) (3.2) 

where the Jacobian V associated with the change of variables from the set of vortex positions 
to the field ip is defined by: 

where e^^^^j^...^^^ is the n-dimensional fully anti-symmetric tensor and summation over re- 
peated indices is implied. The unsigned density, n{l), does not consider the charge of the 
defect and is given by 

n(l) = 5](5(ri-x,(ti)) 

a 

= 6[m]m)\- (3.4) 

If we have products of such densities at equal-times, ^1 = ^2 = t, we write 

p(l)/^(2) = 5(ri - rs) ^ 5(ri - x„(t)) + p(l)p(2). (3.5) 

a 

We use the tildes in (|3.5| ) to indicate that the product p(l)p(2) contains only terms arising 
from different defects. One may write 

p(l)p(2) = 5(ri - r^Hl) + p(l)p(2). (3.6) 

The equal-time charged or signed defect correlation function is given by the average of ( p.6|) 



a(r,t) = (p(l)p(2)) 

= no{t) 6ir) + Gsir,t). (3.7) 



The first term in (|3.7] ) represents self-correlations and 

no{t) = (n(l)) (3.8) 



is the the average unsigned point defect density at time t |[TT|. The second term in ( p.7|) is 

a(r,t) = (p(l)p(2)) (3.9) 



and measures the signed correlation between different defects. It is easy to see that the 
signed defect correlation function, Gg, can also be decomposed as 

Gg = 2Cyy — 2Cya (3.10) 

where C„„ is the correlation function between like-signed vortices and C^a is the correlation 
function between vortices and antivortices. 

We can also define the equal-time unsigned defect correlation function as 

G„(r,t) = (n(l)n(2)) 

= no{t) 6{r)+Gu{r,t) (3.11) 

with 

G„(r,t) = (n(l)n(2)). (3.12) 

As with the signed quantity, one can write Gu in terms of the vortex-vortex and vortex- 
antivortex correlation functions: 

Gu = 2Cyy + 2Cya- (3.13) 



Inverting ( p.lQ ) and ( |3.13| ) and using ( |3.7] ) and ( |3.11| ) one has 



a„(r,t) = ino(t)(5(r) +d.(r,t) (3.14) 
CUr,t) = ^[Gu{r,t)-Gs{r,t)], (3.15) 



where 



a.(r,t) = ^[Ggir,t)+Gu{r,t)]. (3.16) 

As one would expect, for the vortex-antivortex correlations, there is no 5-function contribu- 
tion from self- correlations. 

We are interested in the correlations between different vortices and antivortices in the 
scaling regime and we can define the following scaling forms for the quantities of interest: 
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- to (3.18) 

cu^) = ^£^^jl^^g,(x) + gjx) (3.19) 

c„(i)si^2ji|j|^ = a„W-e.w. (3.20) 

Both Ct,„(x) and C^ai^) are normalized to approach 1 as x ^ oo. 

Since the topological charge is conserved, one has the very important constraint 



J d'^r Gsir,t) = 0. (3.21) 
Using p.7| ) this conservation law can be written in the form 

d'^r Gs{r,t) = -no{t). (3.22) 



Theory, simulation, and experiment indicate that the scaling results ( p. 171) - (|3.2CI|) and 



no{t) = (3.23) 

where A is a constant, hold. Inserting these results in ( p.22| ) leads to the relation 

J d'^x Gsix) = -J. (3.24) 

A measurement of nQ{t) and a choice for L{t) fixes A, allowing one to check that the sum 
rule ( |3.24|) is satisfied. As shown in our theoretical results satisfy this sum rule exactly. 



The results presented above are rather general. To evaluate Gg explicitly one can use 
the gaussian closure approximation, as was done in p[. The evaluation of in this same 
approximation is technically much more difficult than the calculation of Gs because of the 
absolute value sign in the definition ( |3.4| ) of the unsigned defect density n{l). The purpose 
of this paper is to compute G^ and thereby obtain Ct,„ and C^a- In sections VI and VII 
we carry out this calculation with explicit results for n = d = 2. First, however, we must 
briefiy review the calculations of the order-parameter scaling function and the signed defect 
correlation function within the gaussian closure approximation. 



IV. THE GAUSSIAN CLOSURE APPROXIMATION 



Substantial progress has been made in determining the order-parameter scahng function 



T using the theory developed in |12]. In this and related theories, one expresses the order 
parameter '?/!'(r, t) as a local non-linear function of an auxiliary field m(r, t) which is physically 
interpreted as the distance, at time t, from position r to the closest defect. One of the 
physical motivations for introducing m(r, t) is that it is smoother than the order-parameter 
field. Sharp interfaces or well-defined defects produce a non-analytic structure in the order- 
parameter scaling function jF(x) at small-scaled distances x which is responsible for the 



Porod's law decay seen scattering experiments |13[- The expectation, however, is that the 
auxiliary field correlation function f{x), defined below, will be analytic in this same distance 
range. In the case of a scalar order-parameter these expectations are supported by the theory 
||12|| . However, for the simplest theory [Q with n > 1 this is not the case. One finds a weak 
non-analytic component in / and, more significantly, for n = 2 one can trace this non- 
analytic component to an unphysical divergence in the scaling form for the signed vortex 
correlation function Qs{x) at small x. In previous work for n = 2 we showed how these 
problems can be resolved by taking seriously the assumption that the correlations of the 
auxiliary field are indeed smoother than those of the order parameter. We find that it is 
possible to rearrange the theory such that / is analytic in x if we extend the theory to 
include fluctuations about the ordering field and treat the separation between the ordering 
field and the fluctuation field carefully. 

More specifically we can decompose the order parameter ip as 

xjj = a[m] + u. (4.1) 

a is chosen to reflect the defect structure in the problem and is of 0{1). u represents 
fluctuations about the ordering field a and is of at late times. The defect structure 

is incorporated by demanding that a satisfy the Euler-Lagrange equation for the order- 
parameter around a static defect in equilibrium. This determines a as a function of fh. 
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Since we expect only the lowest-energy defects, having unit topological charge, will survive 
to late-times we obtain 

a[rn\ = A{m)m. (4-2) 

In 1^ it was shown that A increases linearly from zero near the defect core and relaxes 
algebraically to its ordered value A = ipo for large m = \m\. 

Evolution under (|2l] ) causes ip to order and assume a distribution that is far from 
gaussian. However, it is reasonable to assume that the the probability distribution for the 



auxiliary field m will be near a gaussian. Indeed, a simple and successful assumption ||I4 
to make is that the probability distribution for m is gaussian with the correlation function 
Co (12) explicitly defined through 

{miil)mji2)) = Sij Coin). (4.3) 

The system is assumed to be statistically isotropic and homogeneous so Co(12) is invariant 
under interchange of its spatial indices. For future reference we also define the one-point 
correlation function 

Soil) = Coin) (4.4) 

and the normalized correlation function 

/(12) = (4.5) 

Since m measures the characteristic distance between defects it is expected to grow as L in 
the late-time scaling regime. This means that Cq and Sq grow as at late times. 

If m is treated as a gaussian variable then its probability distribution is characterized 
by the single function /. This function can be determined by requiring that the equation 
of motion for a be satisfied on average. In |^ it was shown that for n = 2 this require- 
ment produces the following late-time scaling equation for the equal-time order-parameter 
correlations: 



TT 



TT 



X 



V^T + V^J^ - — ^ + — = 



All 2ii 

where, for general ra, JF ( p.3|) is related to / via 

"1 n+ 1 



(4.6) 



n 

2^^ 



1 1 _ n + 2 _ 2 

2' 2' 2 , 



(4.7) 



Here B is the beta function and F is the hypergeometric function [p|,p!5|. In the derivation 
of ( [4.6|) we have defined the scaling length 



L\t) 



2/1 



At. 



(4.8) 



With the theory in this form there will be no leading small-x non-analyticities in the nor- 
malized auxiliary field correlation function / and the small-x divergence in the scaling form 
for the signed vortex correlation function found in earlier theories does not appear. The 
calculation of the scaling form for JF (and thus /) reduces to the solution of the non-linear 



eigenvalue problem ( [4.6|) with the eigenvalue /x. The eigenvalue is selected by numerically 0] 
finding the solution of ( |4.6| ) which satisfies the analytically determined boundary behaviour 
at both large and small x. 



V. THE SIGNED DEFECT CORRELATION FUNCTION 

The calculation of Gs{r,t) (|3.71 ), carried out in 0, begins with the observation that ip 
and m share the same zeros, and that near these zeros we can use (14. 21) to write 



= aom + bom^m + ... (5.1) 

where Oq and Bq are constants which depend on the potential. It is then easy to see that in 
equations ( |3.2| ) and ( |3.4D for p(l) and n(l) we can replace ip{l) with m(l) and the factors 
of Qq and Bq all cancel. Then, assuming m is a gaussian field, it is straightforward to see 
that ^^^(r,^) (|3.9| ) factors into a product of gaussian averages which can be evaluated using 
standard methods. One then finds that Gs{r,t) indeed has the scaling form ( |3.17| ) with 
Qs{x) given by |l|] 
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h{x) 



X 



n-l 



dh{x) 
dx 



with 



h 

7 



_ir_ 

27r 
1 

7r^' 



and 



n 



5(2) ^ l([Vm]2) = 1 
in this theory. The defect density is given by 

rioit) 



-| n/2 



2vr^o(t) 



(5.2) 



(5.3) 
(5.4) 



(5.5) 



(5.6) 



2"/2r(l +n/2) 

Equation ( ^.6|) leads to the result no ~ L^^ which is just a restatement that there is scaling 
in the problem governed by the single length L{t). 

Since / is determined in the theory for the order-parameter correlation function, the 
function Gs(r, t) is fully determined in the scaling regime. The derivatives in (|5.2|) and ( ^.3|) 
make Qs{x) a rather sensitive function of f{x) for small x. 



VI. CALCULATION OF Gu 



The equal-time unsigned defect correlation function G„(r, t) ( p.ll|) can be evaluated 
through a series of steps. We work with general n. The average that needs to be computed 

is iwm 

G„(r,t) = (5[^(l)]|P(l)|5[^(2)]|P(2)|). (6.1) 

The first step is to realize that one can replace by m in (|6.1| ) and write Gu in terms of 
integrals over the reduced probability distribution G(^i,^2): 

(6.2) 
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where 



"^(0 ~ \^IJ-l,IJ'2,-,IJ.n^Vl,V2,-,Vn^l^i^l^2----^lZ. (6-3) 

and 

6) = {5[m{mm{2)\J{5[{i,)l - V,m.(l)]5[(6);: - V,m.(2)]). (6.4) 

The indices and v range from 1 to 2, unless stated otherwise. The gaussian average defining 
^2) is calculated in Appendix A. It is shown there how one can write ^2) in terms 
of the longitudinal and transverse components of the rotated variables (tj)^ (j = 1 or 2) 
defined by 

{tj)l = H (Qp (6.5) 

where M is an orthogonal matrix. We define the longitudinal piece of tj as {tj)1 = {tj)i and 
write {tj)T to denote the transverse pieces: (tj)'^ with /i > 1. One then obtains: 

G{h, = f G^((fi)L, (^2)^) Grmr. {h^ (6.6) 



where 



27lSr 



Ga(ri)^,(r2)L)= (^)"exp-^[ j:m-2U ihh-ihh] (6.7) 

/ ^ \ n(n—l) ~,2 

GTmT,{t2)T)-(^) exp--:^[ EEW-2/tEE(^i)M]. (6.8) 

\27rSTJ 2St ^=2 ;.=2 

The longitudinal quantities 5"!,, /l and 71, are defined in terms of Co and Sq through 

Sl = - ^{C',r (6.9) 

•JO 

= (6.10) 

/l = ^ (6.11) 

7i = (1 - /!)-^ (6.12) 

The corresponding transverse functions are 
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St = S^^^ (6.13) 

Ct = -— (6.14) 

r 

/t = ^ (6.15) 

7l=(l-/|)^^ (6.16) 

Since the matrix M is orthogonal we have 

ViQ=Vit,) (6.17) 

so we may write 

^^(r,^) = / ]ld{t,);d{t,);\V{t,)\\V{t2)\ [^yGLiih)L,{t2)L) GTmT,{t2)T). 



(6.18) 

Under the change of variables 



{to)L = M {s,)l (6.19) 




and, for /i > 1, 



m = M {s,)l (6.20) 
V It 



G„(r, t) becomes 



G„(r,t) = (^^)"^^^(§!?^7Z^^^^^ WtJl) (6.21) 



where 



Wt,/l)= / n^«A^2);:|I^(Sl)||I^(s2)|eXp-^P0i+(^"'2)i-2^ (.-Ol ■ (?2)l] 

X E E([(^i)rJ' + [(^2)^]' - 2/t (si);(^2);:). (6.22) 

Equation (|6.21|) , with the definition (|6.22|) , is one of the central results of this paper. If one 



took this route in evaluating Gs one would arrive at these same equations, only without the 
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absolute value signs in ( |6.22| ). In the next section we will examine the n = d = 2 case and 



go further to derive an expression for that is convenient for numerical work. 



However, first we examine the general expression ( |6.21| ) for Gu in the small-x limit. 
We show in Appendix B that, in the scaling regime, as a; ^ 0, ^ 1 and, surprisingly, 
— s> — 1. To examine this limit for Gu we make the change of variables 

(^i);: = -^r, + \xi (6.23) 



(^2);: = ^0;: - Jx): (6.24) 



for /i > 1 with 



er = 2(1 - /t) (6.25) 



and, for /i = 1, 



(^i)I = -^0L + ixL (6.26) 



(^2)1 = + \xl (6.27) 



with 

€^ = 2(1 + /^). (6.28) 
Then as and go to zero the arguments of the exponentials in ( |6.22| ) go over to 

[(^i);]' + [(^2);]^ - Vt {siXis.X = Kf + Kf (6.29) 

for /i > 1 and 

(si)i + {S2)l - 2h {s,)l ■ {s2)l = [k? + [XL? (6.30) 
for n = 1. The important point is that, as ei and er go to zero, the Jacobians V transform 



as 
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V{s^) = (6.31) 
V{S2) = (6.32) 



Thus the Jacobians differ only by a sign in this hmit and we have 

\V{si)\\V{s2)\ = -V{s,)V{s2) (6.33) 

as X — ^ 0. Since the scahng form for the unsigned defect correlation function Qs{x) is, up to 
a factor of [no(t)]^, given by (|6.21| ) and (|6.22|) without the absolute value signs we see that 



Qu{x) differs from Qs{x) only by a sign 

Gu = -Qs (6.34) 
as a; — s> 0. The relations (|3.19|) and ( p.20|) then lead to the results 



limCt,„(a;) = (6.35) 

and 

lima,(x) = -2e?,(0). (6.36) 

Thus there is a depletion zone at short-scaled distances for like-signed defects. From previous 
work we know that ^s(O) < 0] so the theory gives a non-zero, positive correlation at short- 
scaled distances for unlike-signed defects. These are general results and depend only on the 
gaussian assumption, used to derive ( |6.21|) , and the particular small-x behaviour of / which, 
as is shown in Appendix C, determines the small-x behaviour of the quantities and fi^. 



VII. TWO-DIMENSIONAL 0(2) MODEL 

We have not yet been able to explicitly evaluate ^m(x) for general n = except for small 
and large x. Here we specialize to n = c? = 2. Much theoretical work has recently been 
done on this case and detailed results are available for the auxiliary field correlation function 
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f{x). An additional motivation for examining this case is that simulation and experimental 
results exist for the vortex-vortex and vortex-antivortex correlation functions. 

If we now specialize ( |6.21| ) and ( |6.22| ) to n = d = 2 then the unsigned vortex correlation 



function is given by 



G„(r, t) = [no(t)]^^7ZS?W(/T, h) (7.1) 



with 



WtJl) = J d^siSs2 J{suS2) exp-^[^, + 4-2fL si ■ S2], (7.2) 



where 



J(si,S2) = J d^xd'^y\suX2 - Si2Xi||s2i?/2 - S22Z/i|exp-^[f^ + - 2/t X ■ f ]. (7.3) 

We have simplified the notation by writing Sj in place of (sj)^ and using x and y in place 
of {si)t and (s2)t respectively. The j^^ component of St is written stj. The quantity J is 
evaluated in Appendix B with the clean result: 

J(si, 5-2) = 27r7^|si||s2| J((?) (7.4) 

where 



J{g) =A^l-g + 4V^tan-i (7.5) 



and 



g = fkh ■ hf. (7.6) 

Notice that g depends only on the angle between si and S2 and not on their magnitudes. 
We can then separate ( [7.2| ) into an integration over magnitudes, followed by an overall 
integration over the angular piece: 

/'27r _ rod 2 

AA(/t, Jl) = {^TiflT / dej{g) / ds^ds2 slslexp --[si + sl- 2gL 3,82] (7.7) 

JO Jo I 
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where we have defined 



Ql = fi Si- S2 (7.8) 

and ■§! ■ S2 = cos 6. Expanding the exponential in ( [7.7| ) as a power series in ql and performing 
the integrations over the magnitudes Si and S2 we obtain 



AfifT, h) = {2n)YT /" deJig) f 2^+^ . (7.9) 



£=0 

The symmetries of g and under 6 ^ 6 — 71 allow us to restrict the region of integration 
in (|7.9|) and we use the definition ( p.l8|) for the scaling form to write the final result for Q^^: 



where 5*^, fi and /t are now functions of the scaling variable x. The integral over 6 and 
sum over i have to be evaluated numerically. 

For large x the functions /, fx and are all small and and ( |7.10| ) simplifies since only 



the first two terms in the series need to be retained to give the essential physical features. 
In this limit the integral over 6 is easily performed and one has 



Qu = i+r + -Afi+m--ur- (7.11) 

Since all of /, /t and fi decay as e~^^/^ for large a; 0], rapidly approaches 1 as x increases. 



VIII. COMPARISON OF THEORY AND SIMULATION 

In our previous work we numerically solved the eigenvalue problem ( [4. 61) and deter- 
mined the function f{x) representing the scaling form for correlations in the auxiliary field. 



With this information we can use ( |7.1CI| ) to determine Qu{,x) since for each x at which we 
know f{x) we can calculate /l, /t and 5*^ and perform the sum over £ and then the inte- 
gration over 9. The sum diverges as — ^ 1 (x — and 0^0) but since the smallest x 
we consider is x = .0001 this is not a real problem - we just have to sum up more terms to 
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achieve a set accuracy. The integration over 6 is straightforwardly accomphshed using an 
open Newton-Cotes algorithm. 

Using the relations ( p.l9| ) and ( p.20[ ) we have calculated the results for C^ix) and C^a(x), 
which are presented in Figures |1] and ^ respectively. As expected, Cvv{x) has a depletion 
zone at small-x, which has a characteristic size x ^ 1 {\r\ ^ L(t)). In contrast, C^ai^) shows 
enhanced correlations in the same range of x. We present the results for Qs and Qu in Figures 
^ and Also shown in Figure |^ is the result of the earlier theory ||^,|^ for Qs which displays 
the divergence resulting from neglecting fluctuations. These four figures display the main 
results of this paper. 

Also shown in these figures are the results of simulation 0. The simulation data for the 
scaling forms for the vortex-vortex correlation function and the vortex-antivortex correlation 
function were taken directly from figures 8 and 9 of 0. Relations (|3.19|) and (|3.2CI|) were 



then used to calculate the simulation results for Qs and Qu- As mentioned in the Introduc- 
tion, experimental data also exist for these quantities. The experimental data track the 
simulation results well, but are more noisy and are therefore not shown here. There is only 
one adjustable parameter in all these fits, which is the (unknown) proportionality coefficient 
between the scaling length L(t) used in the theory and that used in simulations |T^. We use 



this freedom to adjust the horizontal scale of the simulation data in Figure |^ to give the best 
match between theory and simulation at intermediate to large-x. The same scaling factor 
is used to rescale the simulation data in the other three figures. It is amusing to see that in 
Figure |^ this rescaling allows the minima of the simulation data and the theoretical curve 
for the earlier theory p[ to coincide. Figures |l| and § show reasonable agreement between 
theory and simulation except for the short distance behaviour of C^- This discrepancy is 
directly related to the behavior of Qg- Since we know that the theoretical expression for 
Qs satisfies the sum rule ( |3.24| ), the question is whether this sum rule is satisfied by the 



simulation results. We have found that the data presented in P lead to the result 

~aJ d'^x Gsix) =0.85±0m, (8.1) 



18 



which is less than the expected result of unity. We are able to compute this quantity directly 
from the data contained in figures 8 and 9 of ||^, using the fact that in units of the scaling 
length they use in these figures no{t) = so A = 1. The error estimate in is due 
to the uncertainties in reading the data, which is somewhat noisy, from the figures in 0. 
What could account for this breakdown in the sum rule? First, there is the possibility of a 
breakdown in scaling and a violation of the scaling relations (|3.171 ) and (|3.23| ). This seems 



incompatible with the simulation results when viewed as a function of time. A second, more 
likely, possibility is that there are some missing vortex-antivortex pairs in the simulation. 
We speculate that there may be a problem keeping track of annihilating pairs at short-scaled 
distances where they may have a very high relative velocity . This may be the source of 



the short- distance discrepancy in Cya between theory and simulation. 

IX. CONCLUSION 

Our work here has concentrated on the statistical properties of point vortices in phase 
ordering systems. The correlations among like-signed vortices found here meet with our 
expectation that vortices with the same charge repel one another at short-scaled distances 
and that screening of this repulsive interaction causes the correlations to fall rapidly to 
zero at large-scaled distances. The case of correlations between unlike-signed vortices seems 
straightforward from the theoretical point of view. Since these pairs are attractive there is an 
increasing probability of finding pairs on their way to annihilation as one goes to short-scaled 
distances. The simulation results seem at odds with this simple physical interpretation. One 
argument around the monotonic behaviour of the vortex-antivortex correlation function is 
that the annihilating pair is speeding up in the late stage of annihilation and therefore the 
probability of finding the pair separated by a short distance is commensurately decreased. 



In recent calculations |T9[, using fundamentally the same theory, one of us found that a 
mechanism already exists in the theory to produce large vortex velocities. These large 
velocities are inferred from a power-law tail in the vortex velocity probability distribution. 
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Bray pO| has shown that this tail results from scaling arguments applied to the the late 
stages of the vortex-antivortex annihilation process. Thus it appears that this speeding up 
process is included in the present analysis. 

In order to better understand the nature of the correlations between vortex-antivortex 
pairs, it would be instructive to work out the joint probability of having a vortex at position 
r with velocity V2 given that there is a vortex at the origin with velocity vi. This calculation 
is under current investigation. 

One of the remaining unresolved questions is: Where are the missing vortex-antivortex 
pairs in the simulations? Further progress is hindered until this discrepancy is understood. 
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APPENDIX A: 

We derive an expression for the reduced probability distribution 

^(^1,6) = (5[m(l)]5[m(2)]n5[(ei);- V,m.(l)]5[(6);:- V,m.(2)]) (Al) 

appearing in the integral formula ( |6.2| ) for Guij.t). We evaluate this gaussian average for 
equal times ti = ^2 = The (5-functions can be represented as integrals and one has 

=y ^^n^^ n^i.<i2MM) e.v-^ik,uQ, (A2) 

where we have defined 

r(gi, 52, fci, k2) = (exp -i[qj ■ m{j) - {kjY^V ^,m^{j)]) . (A3) 
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In these formulae summation over the repeated indices /i, v and j is imphed. The summation 
over j is from 1 to 2, while the summation over and v is from 1 to n, unless stated otherwise. 
Expression (^3) is of the standard form for gaussian integrals 



(exp y dm(\) ■ mil)) = exp ^ y dU2H{\) ■ H{2)Co{12), (A4) 
so a straightforward calculation yields 

+2 n {{ki)X2r, - (^2);gi%) + {c^ - ^) {hr,{k2):r,f, + ^ ihr^ihr,] (as) 

where primes indicate differentiation with respect to r. This expression can be clarified if 
we introduce the orthogonal matrices where 

MI^M: = M;m: = 5,, (A6) 

and 

K = r„ (A7) 
and then transform to the new variables (Wj)'^ defined by 

= mkx- (A8) 

We then obtain 

exp -^[A{q,, gl) + Al{{Wi)l, (H>2)l) + At{{W,)t, (W^W) + A^, ^2, (H>2)l)] (A9) 

where {Wj)1 = (Wj)'^ is the longitudinal part of Wj and (Wj)T is a shorthand referring to 
the remaining transverse parts of Wj ( {Wj)'^ with /i > 1 ). With this choice of variables we 
notice that F can be factored into a transverse and a longitudinal piece. We see that for our 
purposes we are not required to be more explicit than ( |A6D and ( [A?] ) in defining the M^. 
The quantities appearing in the exponential in (|A9|) are 
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gs) = ^0 E ^/ + 2 Co gi ■ (AlO) 

i 

= 5(2) E[(W^i)i]' - 2 C'^ {Wi)l ■ {W2)l (All) 

j 

At{{w,)t, {w,)t) = -2^ E E(w^i);:(w^2);: + ^^^^ e Em,);? (ai2) 

^ /^=2 1^ 11=2 i/j 

AM, q2, {m)L, (H>2)l) = 2 [gi ■ {W2)l - ■ {Wi)l]. (A13) 



One can integrate ( [A^ ) over gl and g2 to obtain 



0/ \" 1 

^ ^ exp--[A'^{{W,)L,{W2)L) + At{{W,)t,{W2)t)] (AM) 



27r5oy 2' 
where we define 

^l((W^i)l, {W2)l) = Sl Y.mj)L? + 2 Cl {Wi)l ■ (H>2)l (A15) 
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with Sl and Cl given by ( |6.9| ) and ( |6.10| ) respectively. Since the transformation (|A8 
orthogonal we have the simple result for G 

G{hM) = GL{{h)L,{h)L) GT((tl)T, (t2)T) (Al6) 



v27r5o/ 

depending on the rotated variable 

{t,)l = M${i,Y^. (A17) 

The longitudinal and transverse parts of t are defined in analogy to those of W . The 
functions Gl and Gt appearing in ( [A16| ) are explicitly given in ( |6.7| ) and (|6.8D. 



APPENDIX B: 

In this Appendix we compute the integral 

J(si,S2) = j d^xd'^y\siiX2 - Si2Xi\\s2iy2 - S22yi\(iW + f - '^fr x-y] (Bl) 

which is needed to evaluate Qui^^it) ioi n = d = 2. The j*'^ component of Si is written Sij. 
To rid ourselves of the absolute values appearing in (|Bl| ) we make use of the identity 
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\x\= r^^(---f]e-^'^'/^ (B2) 



and write 



T/^ ^\ f~^°° dzi dz2 1 d'^ /■ ,2 ,2 1 ^/ ^ ^ ^ .-^^x 

J[si,S2)= — =— = — — — / d xd ?/exp--v4(2;i,Z2,Si,S2,x, (B3) 

J-oo \/ LTl y Iti Z\Z2UZ\aZ2 J ^ 

with 

^2, Si, S2, X, ^) = z\{Sx\X2 - Sx2X\f + ^2("S2iy2 " S22V\f + + - 2/t f ■ ^. (B4) 

The integrations over x and y in ( p3|) are gaussian and so can be readily, if somewhat 
tediously, performed. One has 

= p dz^dz^J (2^)^ .35) 

J-oo V2^V2^ZiZ2 dz^dz2 JD{zi, Z2, Si, S2) 



where the determinant D is given by 

D = il+ zls2')il + zfs,^) - /|(2 + z%' + zls2^ + zlzl[s^ ■ S2]') + f^. (B6) 

The next step is to evaluate the derivatives with respect to zi and Z2 appearing in (|B5D . 
This can be done straightforwardly, and one notices that a change of variables allows one to 
write the integral in ( P5|) as a product of the amplitudes of Si and S2 and an integral whose 
only dependence on Si and S2 is through the dot-product form 

9 = friSi ■ hf. (B7) 

Explicitly, one has 

JisuS2) = 27r^^\s,\\s2\Jig) (B8) 

with 

r^,l_ 7-, |l + 2g+(l-g)(y? + y2-) + (l-9)'CTll ,„m 

J(g) - dy,dy, [l + + ,| + (i _ ' (B9) 

This seemingly complex integral for J{g) is actually pleasingly simple and after some ma- 
nipulations one has 



J{g) = A^l-g + 4v/^tan-i (BIO) 
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APPENDIX C: 



To examine the small-x behaviour of the defect correlation functions we need to know 
the behaviour of /, fx and fi for small x. In the theory we use here, which is discussed in 
detail in H], / is analytic at short-scaled distances and we have 



(CI) 



where a and /3 are constants determined in the theory. To evaluate /t we write Ct (|6.14|) 
in terms of /, function of the scaled length x: 



C = - — — 



1 /' 



X 2da X 

where we have used the result Sq/L"^ = \/2da 0. The result 
definition ( |6.13|) for St and the expansion (|C1|) lead to 



fT = — = l X + ... 



(C2) 

together with the 
(C3) 



for small x. We now examine in the scaling regime and use (|6.9| ) and (|6.1CI|) to write 

1 1 



2/ tl\2 



d 2da 



(C4) 



and 



Cr 



2da 



[f" + fY{f'n 



(C5) 



For small x we again use ( |C1|) and obtain 



h^^ = -l + Oix^). 

The key results here are /t 1 while fL—>-—l when x 
responsible for the relation 



between the signed and unsigned defect correlation functions at x = 



(C6) 



0. This minus sign is ultimately 



(C7) 
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FIGURES 
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FIG. 2. The scaling form Cya{x) for the vortex-antivortex correlation function. The solid curve 
is the result for the theory presented here. The dots represent the simulation data 
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X 

FIG. 3. The scaling form Qs{x) for the signed vortex correlation function. At x = 1 the lower 
solid curve is the prediction of the original theory that does not treat fluctuations, and the 
upper solid curve is the result for the theory presented here which does include fluctuations. The 
dots represent the simulation data [0]. 
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FIG. 4. The scaling form Qu{x) for the unsigned vortex correlation function. The solid curve 
is the result for the theory presented here. The dots represent the simulation data [0]. 
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